Exactly Solvable Fermion Chain Describing a v = 1/3 Fractional Quantum Hall State 
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We introduce an exactly solvable fermion chain that describes a, u — 1/3 fractional quantum Hall 
(FQH) state beyond the thin-torus limit. The ground state of our model is shown to be unique for 
each center of mass sector, and it has a matrix product representation that enables us to exactly 
calculate order parameters, correlation functions, and entanglement spectra. The ground state of 
our model shows striking similarities with the BCS wave functions and quantum spin-1 chains. Using 
the variational method with matrix product ansatz, we analytically calculate excitation gaps and 
vanishing of the compressibility expected in the FQH state. We also show that the above results 
can be related to a v — 1/2 bosonic FQH state. 

PACS numbers: 71.10.Pm, 75.10.Kt, 73.43.Cd 



Introduction. — The fractional quantum Hall (FQH) 
effect is one of the fascinating phenomena in condensed 
matter physics: In a 2D electron system in a magnetic 
field, the Hall conductivity is quantized as o~h — (e 2 /h)v 
with the filling factor given by a rational number v = p/q, 
due to strong electron-electron interactions [Ji- Al- 
though three decades have past since its discovery, the 
importance of this research field is still increasing, partly 
due to new possible realizations of FQH phenomena in- 
cluding flat band Chern insulators |3| and bosonic sys- 
tems of trapped atoms [J] . 

In recent years, there have been theoretical efforts to 
study FQH states in torus boundary conditions which 
can reduce the 2D continuum system in a magnetic field 
to a ID lattice model 0, ||. This approach sheds new 
light on the FQH physics, and is also used to analyze 
new type of FQH states in flat band Chern insulators [7| . 

In this Letter, based on the ID approach, we intro- 
duce a minimal model with an exact ground state which 
describes a v — 1/3 FQH state [Eq. ((3]) below]. We con- 
struct the Hamiltonian in terms of local positive oper- 
ators much like the Affleck, Kennedy, Lieb, and Tasaki 
model for a quantum spin-1 chain |8j. We discuss the 
properties of this model by obtaining exact expressions 
for various correlation functions, order parameters, and 
entanglement spectra. Moreover, excitation gaps are ac- 
curately obtained via variational calculations. 

ID description of FQH states. — We consider 2D inter- 
acting electrons in a magnetic field B on toroidal bound- 
ary conditions, where Li (i = 1, 2) are the circumfer- 
ences of the torus of the corresponding coordinates Xi, 
and Ib = \JhjeB is the magnetic length which will be 
set to unity. As discussed in preceeding works [1,0], the 
system with two-body interaction assumes the following 
ID discretized model, 
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where c\ (cj) creates (destroys) a fermion at site i, and 
the number of lattice sites is given by N s = L1L2/271". 
The matrix-element Vk m specifies the amplitude for a 
process where two particles with separation k + m hop m 
steps to opposite directions. This process conserves the 
center-of-mass coordinate K\ = X)i=i (mod N s ) with 
hi = c\ci, which corresponds to the momentum along x\ 
direction. Therefore, the system with v = p/q can be 
divided into q independent subsystems. 

For the pseudo potential Q which has the v = 1/3 
Laughlin wave function Q as an exact ground state, the 
matrix elements for large L2 are 



^ m oc(fc 2 -m 2 )e- 2 ^ fc2 + m2 ^. 



(2) 



Thus the hopping terms (m ^ 0) are suppressed exponen- 
tially compared to the electrostatic terms (m = 0) in the 
thin-torus (TT) limit L\ — » 0, and the system becomes 
a charge-density-wave state |\)/o) = 1 100100100 ■■ •). In 
order to describe systems with finite L%, we include also 
the leading hopping terms. This expansion in e" ' 1 is 
well controlled, and we expect it to ca ptu re the physics 
also for more general interactions (El EI lid- \2\. 



Model with exact ground state. — Based on the above 
framework, we truncate the long-range interactions of the 
ID model (Q} at v = 1/3 up to the third neighbor (fc + 
m < 3) assuming only \/VioV3o = V21 which is satisfied 
in Eq. fl2j, Via. Then we have 



i=l 



+ otai{c\c\ +z c i+2 c i+1 + H.c.)], 



(3) 



where we have generalized the parameters c%i , , 7$ G R 
to have site dependence. This truncation is valid as an 
expansion beyond the TT limit [l3| . Now we rewrite this 
model in the following form 



n 



[QIQ 



(4) 



2 



where 



OtiCi+lCi+2 + JiCid + 3, 



+ 2- 



(5) 



Eq. (|4]) is clearly a sum of positive operators, thus, the 
spectrum is positive semidefinite (%) > 0. As we will 
show, the (unnormalized) ansatz, 

i*i/ 3 ) = - i*o) = n e ~ u(jt 1*0) . ( 6 ) 



where U = 7»/aj and t/j = i 
unique zero-energy solutions 



provides the 



Note that [U ll U j ] = for 
\i-j\ ^ 2 and LW l+ 2 |*o) = 0. In Eq. ©, the original 
state |- • • 1001 • • •) in l^o) and its "squeezed" counterpart 
| • • • 0110 • • •) cancel each other by acting Qi, and there are 
no next- nearest pairs | - • • 101 • • •). Hence this state satis- 
fies Qi ^1/3) = Pi 1*1/3) = 0, V«, and it is a zero energy 
ground state (Hf ■ Due to the conservation of the center 
of mass, this state has threefold degeneracy for periodic 
systems, even when the parameters have site dependence. 
Our wave function (|6|) gives exact ground states for open 
boundary systems. We can obtain many new zero-energy 
eigenstates at a lower filling than v = 1 /3 by inserting an 
extra anywhere in the root state, |^ ) m @> because 
the insertion of is equivalent to make open boundaries. 
Moreover, if • • • 000101 type configurations are located 
at these "edges", they also give eigenstates with finite 
energies due to the (3f term. 

The uniqueness of the ground state for each center-of- 
mass sector in the v = 1/3 periodic case can be shown 
using the Perron- Frobenius theorem: First, one can show 
that all the states generated by acting with the Hamilto- 
nian on |\l/o) can be reached by successive applications of 
Ui (for different i) [Hj]; thus, the Hamiltonian takes the 
form of a connected matrix in this subspace. Next, with 
a unitary transformation that changes the signs of 7$, all 
the off-diagonal matrix elements can be made negative. 
Now the Perron- Frobenius theorem implies that there are 
no other zero-energy states than (j6]) in this subspace, 
since all its expansion coefficients have the same sign. 
Finally, it follows that all states that are not connected 
to l^o) ( or translations thereof) by the Hamiltonian al- 
ways include finite amplitudes of next-nearest neighbor 
particle which costs energy when 0f > 0; thus, all such 
states have a finite energy. This concludes the proof that 
([6]) is the unique ground state up to translations. 

Correlation functions and order parameters. — From 
the exact solution (|6]) it is possible to calculate rather 
generic quantities such as correlation functions and en- 
tanglement properties of the ground state. For this pur- 
pose, it is convenient to introduce a matrix product (MP) 
representation 16|, |17[ of the ground state wave function 
([6]). In a periodic system with N s = 37V sites, the nor- 
malized ground state wave function ([6]) can be written 
as 



(7) 




FIG. 1: (a) Density functions (hi) in three-sites unit cell, 
and (b) String order parameter Of tring and dual string order 
parameter Og tr j ng as functions of t. OJ tring (OJ tring ) dominant 
in the large (small) t regime plays a role to characterize the 
"superconducting" ("normal") component in an analogy of 
the BCS theory. 



where gj is identified as the following 2x2 matrix, 



fJ.i 



'hi 







(8) 



Here we have introduced the spin-1 representation for 
three-sites unit cell [11]: |010) -> |o), |001) -> |+), and 
1 100) — !• |— ). For the open boundary system, we should 
extract only (1, 1) component of ([7]) instead of taking the 
trace. 

Using the MP formalism for uniform tj = t and infinite 
systems N — ¥ 00, we obtain the density function which 
has three-site periodicity as [see Fig. UKa)], 
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(9) 



This result shows that the density function becomes uni- 
form at t = ±y/2 19]. The single particle correlation 
function is given by (c]cj) — <5y (hi) due to the center- 
of-mass conservation. 

We can also obtain density-density correlation func- 
tions in a similar way. In the infinite-size limit we find 
exponentially decaying correlations 



{ninj)-{ni) 



1 - V4t 2 + 1 



l<-il/3 



1 + V4i 2 + 1 



= e-\- 



-j\/e 



(10) 



As \t\ 00, the correlation length £ diverges which re- 
flects the state |H 1 ) + | — I 1 ). 

Using the suggestive analogy of quantum spin-1 chains, 
we consider nonlocal order in terms of the string order 
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S?) [20(. This is 



string ^ 

known to characterize the Haldane-gap (including Neel) 
state. On the other hand, a "dual" string order param- 



eter 0" trino . = ( e 17r ^=*+i s ?) [2l| can be introduced to 
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FIG. 2: Entanglement spectra {£;}, and entanglement en- 
tropy Sa, as functions of t for (a) infinite-size system and (b) 
finite-size (TV = 2L — 32) system. 



characterize the large-D phases. Using the MP formal- 
ism we find 



(l z = 

^string 



(V4t 2 + 1 - l) 2 
it 2 + 1 
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where TV, j — i — > oo is assumed, and the x, y com- 
ponents are vanishing. In conventional quantum spin 
chains, these two order parameters are usually not finite 
simultaneously. However, the present spin-mapped sys- 
tem breaks the space-inversion and spin-reversal symme- 
tries (e.g., a configuration | 1 — o • • •) occur in (0 but 

| - • • o — h • • •) does not), and also SU(2) symmetry, which 
enables the two different orders to coexist. This results 
consistent with the numerical analysis which concludes 
that the spin mapped v = 1/3 FQH state is adiabati- 
cally connected both from "Haldane" (Neel) and large- D 



phases without closing the energy gap [13, [21 
two string order parameters behave as Oaring 



22] 
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1 for 



and O 



string 



1 for t — > as shown in Fig.[TJb). 
This can be interpreted as the two limits characterize 
"superconducting" and "normal" states in analogy with 
the BCS wave function which is very similar to Eq. ([6]) — 
it has the form of bosonic operators acting on a vacuum 
state. 

Entanglement spectrum and entropy. — We can also de- 
rive the entanglement spectrum (ES) {^} [13] of the sys- 
tem in the spin-1 MP basis via the Schmidt decomposi- 
tion dividing the system into two parts, {k\, . . . , k^} G A 
and {fci+i, . . . , /cat} € B, as 

l*i/ 3 )=AA- 1 / 2 J2 {[9i'--9L] hh [9L + x---9N] hh } 

ii.i2=i,2 

= $>-«*/ 2 |^>^f) (12) 
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FIG. 3: The lowest excitation gaps of the model © at a = 1, 
t — 1/3. The neutral gaps for AK = 1,2 are obtained by 
exact diagonalization (ED) for N s — 27 (dots) and the vari- 
ational calculation with the MP ansatz for N s — > oo (lines). 
E(N S — 1) for 7V S — > oo is also obtained by the MP ansatz. 
Here, /3 2 = 0.77 corresponds to a FQH system on a torus with 
circumference L\ — 5.99. 



where \ipf ) with i = (ji, j'2) are orthogonal states de- 
scribing subsystem A{B), and N^fJ are their norms. In 
the infinite-size limit, L = TV/2 — > 00, one finds 6j6 = 
log(4 + 1~ 2 ) ± log and 6=6 = log(4 + 1~ 2 ) 

[see Fig. EJa)]. The structure of the ES is different from 
that of usual Haldane-gap systems characterized by two- 
fold degeneracy in all ES (6 = 6 and 6 = 6) [? ]■ This 
is because our "Haldane" state is rather close to a Neel 
state which does not have edge spins, due to the lacking 
of the symmetries. We also get the von Neumann entan- 
glement entropy Sa = ^2iCi e ~^ which approaches log 4 
for large t. The finite entanglement entropy is a generic 
property of ID gapped states [24 1. 

In finite systems the entanglement properties are some- 
what altered as shown in Fig. [2jb). Once the correlation 
length £ is of the order of the distance between the cuts 
L, 



the above structure of the ES breaks down |25[ and 
for large enough t the entanglement entropy is instead 
approaching log 2 due to the states |H 1 — • • •) and 

!- + - + •••)• 

Compressibility and excitation gaps. — When we shrink 
the system size as TV S — > TV, — 1 by removing from the 
root state, a 10-type domain wall appears that carries a 
fractional charge [2J e* = e/3 (The fractional charge fol- 
lows from noting that creating three such domain walls 
101010 amounts to adding one electron to the root state 
100100). This excitation energy E(N S — 1) can be ana- 
lytically calculated within a variational approach based 
on the MP formalism. Considering a subspace given by 
© where |* ) is replaced by \^) = |10010100100 • ■ •), 
and appropriate local deformation of the MP state, we 
get a finite value of E(N S — 1) for TV, — > 00 as shown 
in Fig. O Since E(N S + 1) = E(N S ) = as already 
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discussed, we obtain divergence of the inverse compress- 
ibility as expected for the FQH state, 
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The excitation energies in the charge neutral sector 
can also be calculated similarly. We specify these states 
by AK which means the center-of-mass coordinate rel- 
ative to the ground state. Considering AK = 1 and 
AK — 2 subspaces given by ([6]) where l^o) is re- 
placed by states = 1 100100010100 ■ ■ •) and |* 2 ) = 
1 100010010100 ■ ■ •), we get excellent agreement with the 
numerical results of the exact diagonalization which are 
very insensitive to the system size, as shown in Fig. El For 
t = 1/3, the AK = 1 state is the lowest excitation only in 
the very small (3 region, while the AK = 2 state becomes 
the lowest as (3 is increased. The tiny deviation at /3 = 
is due to phase separated states, 1 1010 ••• 0000), which 
have zero energy in this limit. The lowest AK — exci- 
tation is significantly higher in energy. Since the above 
features are qualitatively unchanged from small to suffi- 
ciently large t regions (|t| ~ 1), we identify the AK = 2 
excitation gap as the neutral energy gap of a v = 1/3 
FQH state with toroidal boundary conditions. This re- 
sult is consistent with a recent analysis using the spher- 
ical geometry and the Jack polynomials which identifies 
the neutral gap in the L — 2 angular momentum sector 
26]. 

Bosonic systems. — The present exact argument can 
also be applied to bosonic systems. The Hamiltonian ^ 
with the operators Qi — onbibi + 7i6i_i6j + i and Pi = 
f3ibibi + i defines a v = 1/2 bosonic FQH state that has 
the following exact two-fold degenerate ground state: 
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where \^) = (0101010 •••) and t t = 7 ,/a 4 . In this 
model the spin-1 mapping is also possible as 1 10) — > |o), 
1 02) |+), and |00) |-). Hence, we obtain the 
same effective spin-1 representation ([TJ as in the v = 1/3 
fermion system. 

Conclusion. — We have introduced a ID interacting 
fermion model with an exact ground state that describes 
a v = 1/3 FQH state. We have demonstrated the unique- 
ness of the ground state with periodic boundary condi- 
tions for each center-of-mass sector. We have introduced 
a MP representation of the ground states and obtained 
exact expressions for various correlation functions, order 
parameters and the entanglement spectra. Moreover, the 
excitation gaps have been accurately obtained via varia- 
tional calculations with MP ansatz. We have also shown 
that the present argument can be applied to a v = 1/2 
bosonic FQH state. Although our thin-torus approach 
does not describe a genuine liquid state [llj . it captures 
other important aspects of the FQH physics. 
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